Relaxation of classical many-body hamiltonians in one dimension 



Stefano Lepri 

Max-Planck-Instttut fur Physik komplexer Systeme, Nothnitzer Strafle 38, D-01187 Dresden, Germany 

lepri@mpipks-dresden.mpg.de 

The relaxation of Fourier modes of hamiltonian chains close to equilibrium is studied in 
the framework of a simple mode-coupling theory. Explicit estimates of the dependence of 
relevant time scales on the energy density (or temperature) and on the wavenumber of the 
initial excitation are given. They are in agreement with previous numerical findings on the 
approach to equilibrium and turn out to be also useful in the qualitative interpretation of 
them. The theory is compared with molecular dynamics results in the case of the quartic 
Fermi-Pasta-Ulam potential. 

Keywords: Relaxation to equilibrium, Fourier modes, Fermi-Pasta-Ulam model 
PACS numbers: 05.70.Ln, 05.45.+b, 63.10.+a 
To appear in Phys. Rev. E (1998) 



I. INTRODUCTION AND MOTIVATIONS 

Consider the conceptual experiment where one of the normal modes of an idealized, one-dimensional crystal is 
excited by means of some external source, in such a way that the system is brought far from thermal equilibrium. 
After switching off the external perturbation, it will relax again to the equipartition state, described by the 
canonical measure. The classical question is: how long will it take? Such a basic issue was seriously reconsidered 
in the light of recent results of contemporary nonlinear dynamics, in particular after the discovery that weakly 
nonlinear systems may display extremely long relaxation times. The latter are related to an effective "freezing" 
of some degrees of freedom, i.e. to the slow diffusion in phase space (see for example Ref. Q] for a recent critical 
discussion and further bibliographical references). Because of their simplicity, hamiltonian chains of oscillators 
are suitable model systems to discuss the problem, both from the analytical than from the simulation point of 
view. 

Although some theoretical work indicates that equipartition in a strict sense is always attained in the 
thermodynamic limit for chains at finite temperature (or energy per particle), no explicit estimates of the time 
scales needed are known. The problem has been attached mainly from the numerical side, by looking at the 
time relaxation of suitable indicators of equipartition among the Fourier modes of the chain [§-0. In particular, 
recent studies [|]]8) focused on the dependence of relaxation times on the energy density or temperature, and 
some empirical scaling laws were found. Nonetheless, no quantitative explanation of the latter is sofar given. The 
need for some analytic clue is even more evident as one consider that the computational limits of simulations can 
be rapidly reached. This is of course a particularly serious limitation at very low energies, when the interaction 
among modes and the resulting dissipative effects are extremely weak. 

A further motivation comes from the closely related problem of energy transport in such systems. As is well 
known, the relaxation of fluctuations is strictly connected to transport coefficients, and the existence of slow 
time scales must be reflected somehow in their properties. Indeed, extensive molecular dynamics studies Q 
gave evidence of the divergence of the thermal conductivity associated to the algebraic decay of the correlation 
function appearing in the corresponding Green-Kubo formula. 

The present work aims to give a contribution to the comprehension of both questions by studying the re- 
laxation of chains with an acoustic spectrum. The models we will refer to are introduced Sec. II. When the 
system is not too far from equilibrium, one can rely on well-established theories, like the perturbative and 
mode-coupling approximations that are reviewed in Sect. III. Their validity will be compared with the out- 
comes of numerical simulations for the specific case of a interatomic potential with quartic nonlinearities (Sect. 
IV). Besides the original motivations, the comparison constitutes a nice test of mode-coupling theory and of its 
typical features i.e. the existence of long time-tails with nontrivial exponents, and the nonanalytic behaviour of 
the spectrum of relaxation times. This latter property is actually a peculiarity of noncquilibrium dynamics in 
one dimension. As it will be hopefully clear, this approach will be extremely useful for the problems at hand, 
as they provide estimates of the relaxation times of the Fourier modes of the chain. Rather surprisingly, a 



comparison with some previous numerical results will show that the latter are largely the main cause of slow 
relaxation to equipartition (Sect. V). 



II. ONE-DIMENSIONAL MANY-BODY HAMILTONIANS 



We consider a chain of N anharmonic oscillators and denote with qi the displacement of the Z-th particle from 
its equilibrium position. The hamiltonian reads as is 
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where the usual Born- Von Karman boundary conditions qi = qi+N are assumed and the potential energy is of 
the form 



V(x) = - + V A {x) 



(2) 



with Va denoting the anharmonic part of it. We will consider homogeneous lattices, so that all the masses are 
set equal to unity and pi — q\. The lattice spacing is also set to unity as well as the harmonic frequency, so 
that all variables in the following are adimensional. This also means that the sound velocity, as defined in the 
harmonic approximation, is equal to one. 

In the following we will always refer to the energy density e (energy per particle) and/or to the corresponding 
temperature ksT = Obviously, for a strongly nonlinear systems, the two quantities are not in general 

strictly proportional. 

The (complex) amplitudes Qk of the Fourier modes are defined through the usual transformation 
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Once the hamiltonian is expressed in these new canonical variables, the equations of motion become 
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with Tk being the interaction force among modes, and we have introduced the usual normal-mode frequencies 

(6) 
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III. ESTIMATE OF THE RELAXATION TIMES 



The formulation of stochastic equations for the dynamics of the relevant variables is rather customary to 
describe the relaxation close to equilibrium |Q . The idea is to describe the effective motion of suitable "slow" 
observables by reducing the level of description. The general strategy involves projection on their subspace, 
and results in linear non-mar kovian equations. Whenever a sharp separation of time scales is possible the 
latter reduce to their markovian limit. The memory term determines the relaxation properties and can be 
estimated self-consistently. Our aim is to apply the above procedure to the present system. In the following 
two subsections we summarize the relevant steps. 
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A. General setting 



Due to the conservation law of total momentum, we expect that in the present case the slow dynamics should 
be associated with the long-wavelength Fourier modes Q k with |fc| <C N/2. Moreover, translational invariance 
implies that each mode is uncorrelated from the others so that we can that consider each mode separately. 
Accordingly let us consider the set of Q k and P k as relevant variables and define the projection operator V 
acting on the scalar observable X as 



fx = = £ 



(XQt) n , {xp* k ) 



(7) 
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The projection of the equation of motion leads to 

Qk = Pk (8) 

P fc = -u>lQ k - [ T fe (t - s)P k (s)ds + R k (9) 
Jo 

where R k = (1 — V)P k is the so-called random force which is related to the memory function by the fluctuation- 
dissipation theorem 

T k (t) = 0(R k (t)R* k (O)) . (10) 
and the renormalized frequencies are given, for a generic hamiltonian like (Q), by 

1 1 f e~P v Wdx 

Obviously in the harmonic limit a — > 0, and the usual bare dispersion relation (||) is recovered. Here and in 



the following we will always deal with bounding potentials so that the integrals in ( |ll| ) are convergent. The 
definition ( |TT|) am ounts then to a renormalization of the sound speed from unity to the temperature dependent 
value v = \/\ + a. 

The main object of study will be the normalized correlation function 

Gk{t) = (3u, 2 k (Qk(t)Ql(0)) , (12) 

which is defined in such a way that Q k (0) = 1 . It satisfies the equation of motion pJ] 

Gk +^lGk = - I T k (t~s)g k (s)ds (13) 
Jo 

Introducing the Laplace transforms G k (z) and T k (z) with the definition 

/>oo 

Tk(z) = / e- izt T k {t)dt (14) 



one has that (with G k (0) = 0) 

r (,\ iz + T k (z) , . 

yk(Z) = — T2 . r i x (15) 

z z - uj k ~ izl k {z) 

If the dissipation is small enough with respect to Q k , the transform G k {z) has two poles close to the real axis in 
the complex plane approximatively given by 

±w k -- lim T k (z), (16) 

(at fixed wavenumber and provided that the limit exists) . This corresponds to both a shift of the renormalized 
frequencies and a small damping j k , given by the imaginary part of ( [l6| ) (i.e. by the real part of T k (uj k )). 
Accordingly, the inverse of the latter defines a characteristic relaxation time of each Fourier mode. 



*One can easily convince hymself of this statement by computing for instance the correlation (Qk(t)Qk>(ty) an d imposing 
that (qi(i)qii(0)) depends only on \l — l'\. 
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B. The mode-coupling approximation 



The above results are more or less formal manipulations: we obviously need to compute explicitely the 
memory kernel and the relaxation rates ^ k by resorting to some approximations. A first conceptual difficulty of 
the projection approach is the fact that R k does not evolve with the full Liouvillean operator associated with 
H |10|| . One generally bypass the problem by simply replacing 

(R k (t)Ruo)) - (fk(tmm (17) 

where the last average is on the full Gibbs measure. In such a way, it is also implicitely assumed that slow 
terms possibly contained in T k are negligible in the thermodynamic limit (see below). A second simplification 
amounts to factorize multiple correlations so that the resulting approximate expression of the memory kernel 
T k {z), together with Eq. ([TJ), constitute a closed system of equations for Q k . The latter has to be solved 
self-consistently. 

Let us focus on the Fermi-Pasta-Ulam (FPU) potential 

V A (x) = ^g 3 x 3 + \g^ (18) 



so that T k = jj® + T ( k i] with 



Fk ] =-93U k -j= ^ u kl w k2 Q kl Q k . 2 (19) 
•^fe 4) = -giUk ^2 u kl u k2 uj k3 Q kl Q k2 Q k3 , (20) 

fcl+/C2+&3 — k 

and where the condition on the indexes of the sum is intended to be modulo N (quasi-momentum conservation) . 
Obviously, this represent a reasonable approximation of a generic anharmonic potential in the limit of small 
anharmonicity. In the approximation where multiple correlations factorize in products, the memory kernel will 
be a sum of two terms (averages with an odd number of Q vanish): 

n(t) = /3[<^ 3) (t)4 3) *(o)) + (^ 4) (*)4 4) *(o))] • (21) 

and one can readily evaluate the two contribution from the cubic and quartic terms respectively 

P(4 3 \t)4 3 >(0)) « c 3 ^l £ G kl (t)G k2 (t) (22) 

ki+k2=k 

^ 4) (^f(0)>«^l J2 O kl (t)G k2 (t)g k3 (t) (23) 

k!+k 2 + k 3 =k 

where we have defined the two constants 

c > - 3 wh • c - - 15 WTa? (24 » 

The numerical factors come from the counting all the possible factorizations. Before going further, notice that 
in the present example we can easily understand the assumption of neglecting slow components in (jT^) . The 

force contains indeed a term proportional to Q k \Q k \ 2 /N, which is clearly as slow as Q k itself. As assumed, 
its weight vanishes like l/N, but for any finite N, it may be then regarded as one of the sources of finite-size 
effects in the numerical simulations. 

Even with all the above simplifications the theory remains too complicated to be solved. A further assumption 
is then generally required, namely that the sums in Eqs. (|2^) and (^3j) can be replaced with their values at 
k = JTl|] . For the cubic case this amounts to set k\ = — k 2 = k' in Eq. (]22j), yielding 

P{4 3 \t)4 3 >(0)) « C 3 Lot^J2^(t) (25) 
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For the quartic term we can also extend the summation only to the small k terms that are almost-resonating, 
e.g. those for which ki — — k% w 0. This simplification is justified as in the long-time limit only the slowly 
oscillating contributions should be significative. One can convice hymself that this kind of approximation leads 
to the same result ( p5| ) also for the quartic term, with C3 replaced by C4. Finally, in the limit N — > 00, we let 
2irk/N — > q and replace the sums with integrals. The above hypotesis that we can set k = in Eqs. ( p2|) and 

e amounts therefore to say that the memory kernel can be written in the form T(q, z) = v(z) q 2 for q — > 
. As a result, one gets the self-consistency relation for v from the Laplace transform of Eq. (HQ) 



v{z) « cj™dte^ zt J ^G 2 (q,t) 



(26) 



where C = C3 + C4. This is readily solved by dimensional arguments pJj ], and yields 

u{z) (x -£= . (27) 

This last relation implies that T(q, z) ~ z^ 1 / 3 ^ 2 , so that the limit in Eq. ( [Pf ) leads to a non-analytic dependence 
of the relaxation rates for small wavenumber: 

7(9) oc f^J g 5 / 3 . (28) 

Generally speaking, the behaviour of the relaxation rates with the temperature will depend on the specific 
form of the anharmonic potential. Nevertheless, the q— dependence should be the same for all one-dimensional 
models where the theory applies. 

C. Kinetic vs. hydrodynamic relaxation 

The above self-consistent result is expected to hold for strong anharmonicity and on very long ("hydrody- 
namic") time scales. On the other hand, in the usual perturbative limit, the system is basically a set of weakly 
interacting harmonic oscillators with renormalized frequencies. We can then use the perturbation theory in 
our simple mode-coupling scheme. This basically amounts to neglect the dissipation on the r.h.s. of Eq. ( |26| ) 
and approximate Q(q,t) rs co$u)(q)t. In this limit the factorization of correlations becomes exact. Accordingly, 
Eq. (|2^) reduces then to a simplified version of the usual perturbative formula Jl3[ , where only mode-mode 
contributions are taken into account. It then follows that the spectrum of the dissipation rates is proportional 
to a; 2 times the sum of two terms whose magnitude is given by the constants C3 and C4 respectively. In the 
limit of low temperatures and/or weak couplings the latter scale as 



(.939) 2 £ (from cubic term) 



/3(1+q) 5 / 2 

ffl(i+a)V2 ~ (ff4g) 2 £ 2 (from quartic term) (29) 

for small wavenumbers. It is of basic importance to compare this latter ("kinetic") time scale with the one 
determined in the previous subsection. For example, in the case of the quartic FPU model we have from Eqs. 
(||) and H) 

~ <? 1/3 (54£) 2/3 « 1 • (30) 

This implies that, for small enough 946, the initial relaxation stage is dominated by the kinetic time scale, up 
to some (possibly large) crossover time where the self-consistency effects become relevant. In Section V we will 
come back on this issue to show its importance for the approach to equilibrium. 

IV. MOLECULAR DYNAMICS RESULTS 

As many assumptions are required in the theory, it is important to compare it with the outcomes of numerical 
simulations. To this aim, we considered the quartic FPU case (03 = 0). In this case, the only relevant parameter 
is 546 and g^e <§; 1 correspond to the weakly chaotic regime || . 
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The numerical simulations were performed at constant energy by integrating the equations of motion with 
a third order symplectic algorithm |Q. We generally consider the case where the second constant of motion 
Pq = Qo is identically zero (no uniform rotations of the chain) . Equilibrium initial conditions were chosen either 
by assigning random velocities from a Gaussian distribution at the corresponding temperature, or starting from 
equal mode amplitudes with uniformely distributed random phases. The system is then evolved for a certain 
transient time in order to start the measurements from a more generic phase-space point as possible. In 
computing spectra and correlation functions a Fast Fourier Transform routine has been used, and the data are 
usually averaged over an ensemble of several trajectories (typically between 20 and 200) to reduce statistical 
fluctuations. 

Let us first of all comment on the dynamics of Fourier modes. The correlation of the fluctuating force decays 
on a characteristic time that is expected to be much shorter than the typical period of long-wavelength modes. 
Thus, if we neglect memory effects and assume that a single relaxation time dominates, Eq. (^) reduces to its 
Markovian limit (we come back to the discrete case): 

Qk + ikQk + ^lQk = Rk , (31) 

where now the random force is well approximated by a Gaussian white process 

{R k (t)Rt(t'))=jS(t-t') . (32) 



For simplicity, in Eq. (|3ll) we have also neglected the small frequency shift. Eq. (|31|) explains qualitatively 



the numerical results of Ref. |12| where the dynamics of the modes was studied for the quartic FPU case. In 
particular, the slow diffusion of energy observed there is immediately understood as a consequence of the fact that 
"f k /uj k <C 1 for small k. The renormalization of the frequencies, that was proposed on purely phenomenological 
basis is, in the present context, a straightforward consequence of the projection approach. Furthermore, the 
effective sound velocity can now be explicitely computed by the definition ( jll| ) (at least up to correction from 
Eq. ©). 

At high energy (g^e ^S> 1), where relaxation occours on faster time scales, it is relatively easier to perform 
direct tests of the goodness of Eq. (|3l|). We first verified that the distributions of the real and imaginary parts 
of Qk and Pk are Gaussian. A typical correlation function Gk is reported in Fig. [l| Moreover, we checked that 
the distribution of amplitudes and phase jumps agrees with what predicted from the approximation Eq.(j3l|) 
[ p^| . For completeness, we also measured the dependence of the effective oscillation frequency, and compared it 
with the expected value given by the definition (jfj]) (see Fig. |). 

Clearly, the crucial point to be checked is the energy and wavenumber dependence of the j k . This was 
accomplished by measuring the initial decay of the envelope of Gk (see Fig. |l|), for several es and ks. A very 
good agreement with the mode-coupling prediction is obtained for the dependence of 7^. on the wavenumber. 
The data reported in Fig. [|give a power law with an exponent 1.64, remarkably close to the expected value 5/3. 
Furthermore, Fig. || shows that also the scaling with energy is reasonably obeyed, at least within the limit of our 
simulations. Obviously, the computations become more and more time-consuming with decreasing temperature 
due to the rapid increase of relaxation times and finite-size effects (see below) . This imposes severe constraints 
on the accessible lattice lengths and times. 

A further analysis has been performed on the quantity 

£=Y,E k ; E k = + \co 2 k \Q k \ 2 (33) 

k z A 

which is related to the typical indicators used in equipartition studies. For \k\ <C N/2 (small noise amplitudes), 
we expect from ([H]) that the effective mode energy E k is ruled by the linear equation |l5| 

E k + Jk (E k - (E k )) = R' k (34) 

(see the inset of Fig. Q). This means that the asymptotic behaviour of fluctuations 5£ of the above defined 
quantity are given (in the thermodynamic limit) by 

(6£(t)S8(0)) oc lV e-™* « / * e -r(«)* « h~ 1/2 iort«t c (35) 
V W V " N J 2tt \ i- 3 / 5 for t > t C 



(i 



Here tc is the characteristic time scale at which a crossover between the two relaxation behaviours occours. 
Although its magnitude remains undetermined, we expect it to be very large for small temperature. The 
crossover should become actually observable in the intermediate-energy region. 

The numerical results are in substantial agreement with those predictions. The simulations at g^e = 0.05 
(see Fig. ||) indicate a power-law divergence in the spectrum of 5£ with an exponent -0.5 in the observed 
domain. The statistical accuracy and the length of the simulation are not sufficient to establish whether the 
apparent saturation for co < 10~ 3 is the beginning of the crossover or simply a finite-size effect. A second series 
of simulations in the intermediate energy range i.e. for g^e = 0.45 is reported in Fig. ^. Despite the strong 
finite-size effects at low frequency, the spectra seem to approach a power-law behaviour with the expected 
mode-coupling exponent -0.4 and the data are compatible with (|35| ) with a value of tc of the order of 10 2 . Since 
a further support to the validity of the mode-coupling results is also reported in p6| , we reasonably conclude 
that no significative deviations from the theory itself are observed in those type of measurements. 



V. THE PROBLEM OF RELAXATION TO EQUIPARTITION 

At this point we want now to discuss some consequences of the above analysis on the so called FPU-problem. 
Let us first focus again on the case of a purely quartic nonlinearity that has been intensively (re)studied by 
several authors in recent years ^,^,^|. The numerical experiments have been performed by feeding the initial 
energy in a packet of modes and looking at the decay in time of suitable indicators of equipartition (see the 
quoted references for details). 

In the case when the initial excitation is around a long- wavelength mode of wavenumber q,, and the system 
is not too far from equipartition, we expect, from the discussion of Section III, that the slowest time scale will 
be of the order of 



/ (54 eq*) for 34E < 1 /q^\ 

TE = \ I W/4 -5/3 , ^ 1 ■ W 

(54 e) 1/4 (?* for g 4 e > 1 



The subscript E is precisely meant to remark that it refers to the fluctuations close to the equilibrium state. 
The behaviour in the high-energy limit is obtained by means of Eq. ( pq ) taking into account that e is roughly 
proportional to the temperature, as well as the fact that (1 + a) grows with the square root of the temperature 
itself (see Eq. (0)). 

Indeed, the predicted regimes have been numerically observed [0. A convincing numerical evidence of scaling 
laws ( |36| ) has been recently reported in Ref. ||. This seem to suggest that the linear theory already captures 
the quantitative features. Moreover, qualitatively similar results were found for different potentials in Ref. |Tt|| . 

For what concerns the case of larger , the simulations reported in Ref. also showed that relaxation can 
be one or two orders of magnitude faster than in the previous situation. Even if the theory presented here 
applies better to the case of small q* , it is worth to remark that its consequences are also consistent with this 
observation. 

Obviously, the very initial stage of the approach to equilibrium may well occour on a different time scale. In 
the present framework, the latter should be interpreted as a real "partial equilibration" time scale tne- On the 
other hand, it seem rather reasonable that t^e will be strongly dependent on the chosen class of initial conditions 
and on the specific form of the potential. As a matter of fact, tjvb is determined by a pure nonequilibrium 
dynamics and is of course unaccessible to the linear theory presented here. Despite of this, the present results 
are useful to identify this initial stage. As an example, we are now able to make an instructive comparison with 
the numerical estimates of t^e- In Ref. || it is in fact found that, for long- wavelength excitations (q» <C 1), 
t ne ~ N 1 / 2 {q*e)~ x , so that 



tne 
te 



rA 1 / 2 ^ « 1 (37) 



for small q» and/or small energy. We can then conclude that the equilibration process is, at least for this class 
of initial conditions, mainly dominated by the linear regime. In other words, we can naturally understand 
it as an initial fast relaxation to a quasi-equilibrium state followed by the slow diffusion of energy from the 
long- wavelength modes. 

As already mentioned, the mechanisms determining tne may be however rather complex and of very different 
nature depending on the initial state. An example, is the case of zone-boundary initial conditions rs 7r). 
Those rapidly decay into localized chaotic excitations |L9]|| , whose lifetime mainly determines the time to reach 
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the quasi-thermalized state. Nevertheless, the comparison of our results with the numerics may indicate that 
such lifetime is, even in this case, considerably shorter than that of relaxation of Fourier modes 

The high-energy scaling of Eq. (|3^) can be generalized to an algebraic potential of the form Va [x) = g n x n , 
with n being an even integer. In this case, by extending Eqs. (53) and (Eq), it is found that te is proportional 

to £ (V«-V2). 

Finally let us comment on the the FPU model with purely cubic potential (g± — 0) . A recent numerical study 
IjJI showed a clear evidence for a divergence as e -3 of the relaxation time. Nevertheless only chains as short 
as N — 32 were considered, and it is not completely clear if the systems is above the equipartition threshold 
H . A more detailed analysis would of course be desirable to check the dependence of this scaling on N and to 
compare the results with the ones presented here. 



VI. CONCLUSIONS 



The simple mode-coupling approach provided a valuable amount of qualitative and quantitative information 
on the relaxation times of the Fourier modes in an hamiltonian chain as ([j]). In particular, a clear physical 
interpretation of them can be achieved. For the paradigmatic example of the quartic FPU model, the crossover 
between the two scaling regions at g^e ~ 0.2 (see again Fig. ||), could be regarded as the temperature threshold 
beyond which self-consistency and "hydrodynamic" effects play a major role. Remarkably, such a scale is roughly 
equal to the so-called strong stochasticity threshold (g^e « 0.1) above which new dynamical effects (i.e. faster 
diffusion in phase space) are believed to appear . From this point of view it would be challenging to try to 
connect the present results with strictly dynamical properties. 

The explicit estimates given in the present work turned out to explain several previous results on the approach 
to equilibrium at finite temperature. As esemplified above, they are a precious information for the numerical 
study. For example, the slow temporal decay of the equipartition indicators may be better understood having 
realized the existence of the long-time tails like those of Figs. || and [| Indeed, it is clear that the slow diffusion 
of energy has to be taken into account when studying the relaxation from an arbitrary initial condition. More 
generally, such time scales should be always considered in practice when performing simulations with chains of 
large sizes. 

Although the existence of long-time tails is not surprising for a low-dimensional system, the substantial 
agreement between the theory and numerical simulation is nonetheless a relevant result by itself. Actually, the 
direct verification of mode-coupling theories in one dimension is still current object of study, for example in the 
field of lattice gases |Q. Furthermore, the validity of such theories is not granted in general. It is in fact known 
that they fail in predicting the characteristic relaxation times of spin-waves of the Heisenberg model in one and 
two dimensions pl| . 

Finally, let us remark that the theory allows to estimate the long-time tails of the heat-flux correlation, which 
is directly related to the thermal conductivity of the chain. This allows to explain quantitatively the divergence 
of such transport coefficient observed in the numerical simulations pjlq] . 
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FIG. 2. Comparison between the numerical and expected renormalized frequencies. The symbols are obtained 
measuring the oscillation period of Q\ for N = 64, 128 (diamonds and triangles respectively). The solid line is the 
theoretical value, Eq. ( pd| ) and the dashed one is the empirical formula of Ref. [12]. The small systematic deviations are 
due to the small frequency shifts expected from ([l6|). 




FIG. 3. The wavenumber dependence of the relaxation rates 7(g) at g^e = 8.8 for the quartic FPU potential. All the 
points were obtained from the initial decay of the envelope of Q\ for increasing values of N up to N = 2048. The solid 
line is a power-law fit q 1 ' 64 . 
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FIG. 4. The scaling of the relaxation rates with the energy density e for N = 64, 128, 256 (circles, squares and pluses 
respectively). The inset shows the low energy part. Straight lines correspond to e 2 and e 1//4 respectively. 




FIG. 5. Low-frequency part of the spectrum of S£ for the quartic FPU model for g^e — 0.05. The solid lines refer to 
simulations with N = 2048,4096, the dashed line correspond to a lu^ 1 ^ 2 . A power-law fit gives an exponent -0.49±0.01 
for 10~ 3 < uj < 10~ 2 . 
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